Anisotropic Inverse Cascade toward Zonal Flow in Magnetically Confined Plasmas 
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We propose a new mechanism for the generation of zonal flows in magnetically confined plasmas, 
complementing previous theories based on a modulational instability. We derive a new conservation 
law that operates in the regime of weakly nonlinear dynamics, and show that it serves to focus the 
inverse cascade of turbulent drift wave energy into zonal flows. This mechanism continues to operate 
in the absence of the separation of dynamical scales typically assumed in instability calculations. 
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Zonal flows refer to a class of highly anisotropic flows 
that emerge spontaneously in response to nominally 
rather weakly anisotropic trends in the environment. 
Examples in geophysical fluids include strongly sheared 
east-west jets in planetary atmospheres (Jupiter, in par- 
ticular). North-south variation of the Coriolis parameter 
provides an obvious underlying anisotropy, but is far too 
weak to directly explain the observed flow patterns. 

Similar flows are believed to exist in magnetically 
confined quasi-neutral plasmas, e.g., tokamaks [1], re- 
sponding in this case to gradients in the magnetic field, 
background ion concentration, plasma temperature, etc. 
Zonal flows take on added importance in plasmas because 
they are believed to provide transport barriers, leading to 
the low to high confinement (L-H) transition, and hence 
may aid the goal of controlled fusion. 

There have been many calculations elucidating condi- 
tions under which small-scale drift wave turbulence can 
produce a modulational instability, leading to exponen- 
tial growth of a zonal flow pattern [2|-|6[ . The calculations 
typically assume a large separation of scales, enabling a 
simple description of the growth of an existing zonal flow, 
pumped by sufficiently strong resonant small scale fluc- 
tuations. Broader conditions can probably be formulated 
in terms of the shape of the drift wave spectrum 0. 

The goal here is to show that energy transfer from 
small scale turbulence to large scale zonal flow is a general 
physical phenomenon in plasmas, operating irrespective 
of whether the dominant interactions are local or nonlo- 
cal in scale. We derive a new "extra" conservation law, on 
top of those for energy and momentum/enstrophy. The 
inverse cascade of wave energy follows from standard ar- 
guments involving balance of energy and enstrophy flux 
in the spectral domain 0. The extra invariant places 
further constraints on the energy flux, forcing it more 
and more strongly into the zonal wavevector sector with 
increasing scale. This provides a general mechanism for 
the observed amplification of zonal anisotropy, and will 
be demonstrated quantitatively through analysis of a cor- 
responding spectral function 0k- 

A similar invariant exists in the quasigeostrophic (or 
CHM [|| ) equation [t| [l3| , and in the shallow water sys- 



tem [12j | of geophysical fluid dynamics. However, the 
plasma system is significantly more complicated, produc- 
ing, for example, in addition to the usual CHM "vector" 
nonlinearity, a "scalar" nonlinearity 0, Il3l - [l5| . We find 
it unlikely, for example, that the generalized Hasegawa- 
Mima (GHM) [l3[ equation (with both nonlinearities) 
possesses an extra invariant — this issue will be discussed 
further below. Given the added levels of approximation 
entering such reduced equations, we base our derivation 
directly on the more general effective two-dimensional hy- 
drodynamic equations from which they are derived. We 
account for (smooth, large scale) inhomogeneity in the 
electron temperature (which leads to the scalar nonlin- 
earity), applied magnetic field B, fluid pressure P, and 
background ion density no, all on the same footing. Fur- 
thermore, we do not assume any particular common di- 
rection of variation of these parameters (e.g., the "ra- 



dial" direction in a tokamak [16|). A single direction 7 
emerges naturally, namely the local gradient of the ratio 
no/B, with zonal direction orthogonal to it. 

The effective two-dimensional hydrodynamic equations 
for a magnetically confined plasma, in the xy-plane nor- 
mal to the applied magnetic field B = Bi, are 



<9 t v + (v • V)v + /z x v = -V$ 
d t n + V • (nv) = 0, 



(1) 



where n is the ion number density and v = (v x ,v y ) the 
velocity. The ion cyclotron frequency is / = qiB/rrii, 
where qi,rrii are the ion charge and mass. The potential 
<t> = (qi/rrii)tp + W(n) consists an electric potential ip 
(with E = — V<p), and an ion pressure term (W being 
the solution of rmnW' '(n) = P'(n), with pressure P as- 
sumed a function of the density alone). Electrons move 
rapidly along the magnetic field lines, and are assumed to 
be in local equilibrium at a temperature T. The density 
then follows the Boltzmann distribution, n = noe eip / T 
(so no is the ion concentration when <p = 0). The fluctu- 
ation contribution to B from currents generated by v is 
assumed small compared to B, and is neglected. 

We will see that the slow drift wave modes in {T]) are 
weakly coupled to all other motions, and hence support 
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a separate set of approximate conservation laws. Most 
importantly, the drift wave dispersion law, essentially 
uniquely, admits the new conservation law [llfl . 

Irrespective of the relation between $ and n (if any), 
((T|) leads to convective conservation of potential vorticity 



independent fields, we therefore consider a supplemented 
17] invariant 



(3 t + vV)Q = 0, Q = (C + f)/n, 



(2) 



where £ = V x v = d x v y — d y v x is the vorticity. In 
the weakly nonlinear limit one may approximate Q w 
(Q + f)/n , where Q = £ — fh and h = (n — n n )/n . 
From (JXJ) one obtains the (exact) equation of motion 

d t Q + V- (Qv) = /(l + fc) 7 . v, 7 = Vln(n //). (3) 

Equation ([2]) gives rise to the usual infinite hierarchy 
of conserved integrals. We derive here an invariant of a 
different type. It is quadratic in Q, 

I=\ J d 2 n J dV^n.rajetn^etra,*), (4) 

with some (symmetric) kernel X to be determined. Drift 
wave energy and momentum may be expressed this way, 
but for the plasma system ([1} there is an extra choice. 

We utilize two small parameters. First, we assume 
that no,/, T vary slowly, namely on the scale 1/7, 
which itself is much larger than dominant scale L of 
variations of the fields Q,v. The quantity p = jL 
is the small inhomogeneity parameter. The quantities 
|Vno|/no, |V/|//, |VT|/T, |V 7 |/ 7 are all 0(p/L). This 
is similar to the usual beta-plane approximation in geo- 
physical fluid dynamics, and will allow us, at a critical 
stage, to perform a local Fourier analysis. Second, weak 
nonlincarity constrains the dimensionless characteristic 
field amplitude A ~ Q/f ~ v/fL. The divergence term 
in (J3)) should be much smaller than the linear term on the 
right hand side, leading to the small nonlinearity param- 
eter e = A/ fx. More generally, for a complex turbulent 
state, this parameter should be small on all length scales, 
not just the dominant scale L. 

Using typical D-T fusion plasma parameters, T ~ 
10 KeV, B ~ 5 T one obtains Larmor radius p = 
~ 3 mm. The drift velocity is v 
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(p/ LB, where tp ss Th/e is estimated from the Boltzmann 
relation, yielding (/f ~ h(p/L) 2 , Q/f ~ [1 + (p/L) 2 }^ 
and hence e ~ Q//L7 ~ [1 + (p/L) 2 ]h/p. Using in- 
homogeneity scale I/7 ~ 1 m, density fluctuation scale 
h ~ 10~ 2 , and zonal flow scale L ~ 10 cm, one therefore 
obtains p, e both of order 10 _1 . These are indeed small, 
well within the range of validity of the theory to follow. 

We will prove that there are only three independent 
choices of the kernel X, for which I is approximately 
conserved, i.e., may be considered constant over very long 
time scales (made more precise below). The simplest 
way to do so would be to bound \I\/I. Unfortunately, 
/ contains oscillatory terms that have small amplitude, 
but whose time derivatives do not. Taking v, Q as the 
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d 12 QlF 12 -V 2 + | / <h:2:\Q\ C2-M12:; • V.i 



g / rfl23il23QlQ2fi3, 



(5) 



in which, to condense the notation, numerical subscripts 
stand for the argument: d 12 = d 2 r 1 d 2 r2, v 2 = v(r 2 ,£), 
M 123 = M(rj, r 2 , r 3 ), etc. For small p, e the added terms 
will be found to be much smaller than / but their time 
derivatives generally are not, and the kernels X, F, M, Y 
will be determined by demanding that \I\/fI be small 
[l8j |. The kernel M is symmetric in its first two argu- 
ments, and Y is symmetric in all three. Other cubic 
terms are possible, involving different combinations of 
<2,v, but, due to the structure of ([1]), turn out not to 
contribute, so we drop them from the outset. 

We will see that Fi 2 = O(p), hence to compute 1 it 
suffices to approximate (|T|) , (J3j) by 

d t v = -fz x v - V$, $e = /7-v- V- (Qv). (6) 

Here $ depends only on the local density, and is taken 
to vanish for the steady state plasma. It suffices as well 
to use its linearized form 



* = Th = (T/f)(C-Q) 



(7) 



with the slow function T(r) = T(r)/m, in the standard 
cold ion limit where one neglects the pressure term. 
We define for convenience the combinations 

Gi 2 = (T 2 // 2 )V 2 • F12, #123 - (T 3 // 3 )V3 • M123. (8) 

Using ([T]) and integrating by parts where necessary 
to remove spatial derivatives from the fields, and collect- 
ing terms, one then obtains, to requisite order, I with 
the same four terms as in ([5|), but with corresponding 
(appropriately symmetric) kernels, 

X12 = —{G12 + G21) 

P12 = X 12 hl2 + z x P 12 / 2 + V 2 x G12 

M123 = £ 2 3V 3 ^13 + <Sl 3 V 3 ^ 2 3 + Z X Mi 23 / 3 

+ V 3 x #123 + ^123/373 

^123 = -(#123 + #132 + #321 ), (9) 

where 612 = <H r i — r 2), V x G = (d y G, —d x G), etc. 
A number of terms of higher order in p, e have been 
dropped. These physically represent higher order non- 
linearity, including interactions between drift waves and 
other modes contained in ([1]). The vanishing of Xi 2 and 
^123 produce the antisymmetry conditions 

Gl 2 +G 2 i=0, #i 23 + #i 32 + # 32 i = 0, (10) 
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while the vanishing of F12 and M12 produce 
F12 = 7-V2G12 + (z x 72)^12 

/2 



M 



12.5 



-V 3 ifi 



23 



73^123 
1 



13 



J13V3X 



23 



(11) 



Substituting (fTTT) into the right hand sides of (O, one 
obtains the closed equations 



C2G12 = (72 x V 2 )^12 
^3^123 = (73 X V 3 )Yi23 
<5l3 



V3# ) x (VaJfis 



V 3 — 

/a 



/3 

X (V3X23). 



in which we have defined the operator 

1 



C 



f 



I 
7" 



(12) 



(13) 



We may write formally G12 = fC^Xi^, where K, = 
£ _1 (7 x V). Since K2 depends only on r 2 , it does not 
have any particular symmetry under interchange of 
and r 2 . Therefore, given that X 12 = X 2 i is symmet- 
ric, it is generally impossible to enforce antisymmetry 
of Gi2- However, to leading order, one may ignore the 
r dependence of all the parameters, replacing them by 
constant characteristic values. In this case K, is trans- 
lation invariant, and one seeks translation invariant so- 
lutions X12 — X{ri — ra), with X(r) even. Hence 
G\2 = G(ri — Y2) is odd, and the first of equations (|T0|) 
is automatically satisfied. 

Similarly, given the symmetries of X12 and Y123, it 
is impossible to enforce antisymmetry of K123 beyond 
leading order. However, again replacing all parameters 
by constants, a consistent solution for Y123 does exist. It 
is most conveniently expressed in Fourier space, where, 
in particular fK. — > ifi(k), in which 
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fl x k 

k 2 + a 2 



P 



T 



exhibits the usual drift wave dispersion relation \li 
obtains: 



123°123 



123 



r» 2 + n 3 ) 
1 
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1 



fe 3 



^1 



k 2 



k 3 



a' 



1 



cr 



hi 



X 2 



(14) 
One 

Xi 
(15) 



with <5i23 = (27r) <5(ki+k2+k3) reflecting translation in- 
variance, and A123 = ^ki xk 2 = ^k 2 xk 3 = ik 3 x k 2 be- 
ing the area of the resulting triangle formed by ki, k2, k 3 . 



Up to now, X(k) is arbitrary. However, (fTS"]) displays 
a divergence on the "three-wave resonant" surface 

k 1 + k 2 + k 3 = o, n(k x ) + n(k 2 ) + n(k 3 ) = o. (16) 

Only if the term in square brackets vanishes on this sur- 
face does a nonsingular Yx23 emerge, and this places 
rather stringent conditions on X, that we will now dis- 
cuss. On this surface, the vector /ki2 3 = kif^ — k2^i = 
k2^ 3 — k 3 il2 = k 3 fii — kif2 3 is also symmetric under 
cyclic permutation of its indices. Defining the "zonal" 
and "radial" wavenumbers p = —7 x k, q = 7 • k, one 
may write the term in square brackets as 



7 x ki 



23 



P1P2P3 



pxXx + p 2 X 2 + p 3 X 3 



(17) 



and from its vanishing one therefore obtains the condition 
that pXk also be conserved on the resonance surface. 

The set of kernels satisfying this condition has been 
investigated at length [j| [l(| ■ In addition to the obvious 
choices X^ — 1 (enstrophy/zonal momentum), — 
fl(k)/~ffp = (k 2 + a 2 ) -1 (energy), the extra invariant 



1 

P 



X^f = — arctan ■ 



i(q +pV3) 
k 2 



1 , a(q-py/3) 

- arctan -5 . 

p k z 



(18) 

The existence of the corresponding invariant (j4]) in plas- 
mas is the fundamental result of this paper. Since so- 
lutions to (fT0|) do not exist beyond leading order in (j,, 
neither does extra conservation. Intuitively, higher order 
accuracy in fj, must account for VT, but X M depends 
only on the one direction 7 (via definition of p, q). The 
GHM equation [l3| (with comparable vector and scalar 
nonlinearities) has both gradients, and it appears impos- 
sible to find an extra invariant at all (even if VT || 7) 
20]. Only when one nonlinearity dominates, and cor- 
respondingly, one of the gradients can be disregarded, 
does an extra invariant emerge. Our results show that a 
consequence of ([H is that VT is always higher order. 

To understand the significance of ([TBI , and its rela- 
tion to the formation of zonal flows, consider the linear 
combination 1 = 1 — 2\fiaE = f (jj^E^d 2 k / (2tt) 2 , where 
Ek = X^ |Qk| 2 is the energy spectrum, and 



8V3a 4 p 2 (g 2 +p 2 /5) 
P 



2V3 

0[(a/kf] 
0[( P /q)% 



(19) 



k/a > 1 
p/q < 1. 



Remarkably, X E and X M have identical (up to a factor 
2\^3a) asymptotic behavior for large k and for small p, 
so <j> — > when k — > 00 or p — » 0. The inverse cascade fol- 
lows from the spectral balance needed to maintain conser- 
vation of energy and enstrophy 2l|. As illustrated in Fig. 
[T] the extra conservation provides additional constraints 
leading to zonal flows. Specifically, it follows from (TT9"|) 
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FIG. 1: 3D plot of the ratio (|19[) . which measures how much 
extra invariant I is attached to a unit amount of energy with 
wave vector k. This ratio is even in both "zonal" wavenumber 
p — —7 x k and "radial" wavenumber q = 7 ■ k, with the 
direction 7 defined in Q, and is plotted for 0.4 < k/a < 2. 
If one considers transporting a unit of energy from small scales 
(large k) toward the origin, one is forced to navigate around 
the peak, and squeeze into the valley along the g-axis. 

that a unit of energy at large k/a carries very small <j>. 
Therefore transferring energy towards the origin (either 
by local cascade, roughly along a level curve of 0, or di- 
rectly in a single jump) requires correspondingly small 
values oip/q: it must squeeze around the peak into the 
valley along the y-axis. The resulting flow is zonal with 
velocity along the p-axis, i.e., orthogonal to 7. 

This conclusion is very general, following from a robust 
conservation law [22| that operates under most situa- 
tions considered in the literature. Bounds on spectral en- 
ergy transport imposed by (|18p should help inform future 
more detailed flow computations. Care, as well, should 
be taken in analyzingreduced models obtained from ([1} . 
The CHM equation |8| possesses the extra invariant, as 
well as enstrophy and the infinite potential vorticity hi- 
erarchy. However some generalizations of this equation 
fail to do so, at least in certain parameter ranges. Ad- 
ditional conservation might be restored by re-including 
some neglected terms. 
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